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ON I Results from our recent investigations of quark number susceptibilities in both quenched 
and 2-flavour QCD are presented as a function of valence quark mass and temperature. 
A strong reduction (~40%) is seen in the strange quark susceptibility above T c in both 
cases. A comparison of our isospin susceptibility results with the corresponding weak 
coupling expansion reveals once again the non-perturbative nature of the plasma up to 
3T C . Evidence relating the susceptibility to another non-perturbative phenomena, pionic 
screening lengths, is presented. 

1. INTRODUCTION 

43 . 

Investigations of quark number susceptibilities from first principles can have direct 
Ph! experimental consequences since quark flavours such as electric charge, strangeness or 
baryon number can provide diagnostic tools for the production of flavourless quark-gluon 
plasma in the central region of heavy ion collisions. As emphasized by B. Miiller j3|, 
and Koch in this conference, fluctuations in conserved charges can act as probes of 
quark deconfinement since they are different in the hadronic and plasma phases in simple 
models. Indeed, excess strangeness production has been suggested as a signal of quark- 
gluon plasma almost two decades ago ||. Lattice QCD can provide a very reliable and 
robust estimate for these quantities in both the phases since in thermal equilibrium they 
are related to corresponding susceptibilities by the fluctuation-dissipation theorem: 

{6Q2) k = xq{t ^ q = 0) • (1) 

Here fiQ is the chemical potential for a conserved charge Q, and Z is the partition function 
of strongly interacting matter in volume V at temperature T. Although Lattice QCD is 
unable to handle finite chemical potential satisfactorily at present, and cannot thus yield 
any reliable estimates of a number density, the susceptibility above, i.e., the first derivative 
of the number density at zero chemical potential, can be obtained reasonably well using 
the conventional simulation techniques. 

Quark number susceptibilities also constitute an independent set of observables to probe 
whether quark-gluon plasma is weakly coupled in the temperature regime accessible to 
the current and future planned heavy ion experiments (say, 1 < T/T c < 10). A lot of the 
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phenomenological analysis of the heavy ion collisions data is usually carried out assuming 
a weakly interacting plasma although many lattice QCD results suggest otherwise. It has 
been suggested that resummations of the finite temperature perturbation theory may 
provide a bridge between phenomenology and the lattice QCD by explaining the lattice 
results starting from a few T c . As we will see below, quark number susceptibilities can 
act as a cross-check of the various resummation schemes. Earlier work on susceptibilities 
[0] did not attempt to address this issue and were mostly restricted to temperatures very 
close to T c . Furthermore, the quark mass was chosen there to vary with temperature 
linearly. We improve upon them by holding quark mass fixed in physical units (m/T c = 
constant). We also cover a larger range of temperature up to 3 T c , and the accepted range 
of strange quark mass, in our simulations. 



2. FORMALISM 

After integrating the quarks out, the partition function Z for QCD at finite temperature 
and density is given by 

Z = J VUe~ Sg det M(m u , /J, u ) det M(m d , /j, d ) det M(m s , pi s ). (2) 

Here {U v (x)}, v = 0-3, denote the gauge variables and S g is the gluon action, taken to be 
the standard Wilson action in our simulations. Since we employ staggered fermions, the 
Dirac matrices, M, are of dimensions 3NgN t , with N s (N t ) denoting the number of lattice 
sites in spatial (temporal) direction, rrif and \Xj are quark mass and chemical potential 
(both in lattice units) for flavour /, denoting up(u), down(d), and strange(s) above. 
The chemical potential needs to be introduced on lattice as a function g{p) and g{—\i) 
multiplying the gauge variables in the positive and negative time directions respectively, 
such that M i) <?(//) ■ g(—fi) = 1 and ii) the correct continuum limit is ensured. While 
many such functions g can be constructed, exp(/x) being a popular choice, the results for 
susceptibilities at \i = can easily be shown to be independent of the choice of g even 
for finite lattice spacing a. From the Z in eq. (0), the quark number densities and the 
corresponding susceptibilities are defined as 



TdlnZ dn f T 



d 2 Z 1 dZ 1 dZ 



(3) 



V d\if dfif V [Z djjfdjjf Z dfif Z dfxj 

To lighten the notation, we shall put only one subscript on the diagonal parts of \. 

In order to obtain information for quark-gluon plasma in the central region, we evaluate 
the susceptibilities at the point fif = for all /. In this case, each n./ vanishes, a fact that 
we utilize as a check on our numerical evaluation. Moreover, the product of the single 
derivative terms in eq. (^) vanishes, since each is proportional to a number density. We 
set m u = m d < m s . Noting that staggered quarks have four flavours by default, Nf = 4, 
and defining fi 3 = jj, u — fi d , one finds from eq. (||) that the isotriplet and strangeness 
susceptibilities are given by 

Xs = l~Ox{m u \ Xs = ^Pi{m s ) + \o 2 {m s )} , (4) 

where 1 = (Tr (M"M _1 - M' M '.\/'.U )}. 2 = ((TrM'M -1 ) 2 ), M' = dM/dpt and 
M" = d 2 M/dfi 2 . The angular brackets denote averaging with respect to the Z in eq. @. 
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One can similarly define baryon number and charge susceptibilities. We refer the reader 
for more details on them to Ref. 0. 

In the discussion above, quark mass appears as an argument of O and implicitly in 
the Boltzmann factor of Z. Let us denote it by m va i and m sea respectively. While the 
two should ideally be equal, we evaluated the expressions above in steps of improving 
approximations (and increasing computer costs) by first setting m sea = oo for all flavours 
(quenched approximation |L|) and then simulating two light dynamical flavours, by setting 
m sea /T c = 0.1 (2-flavour QCD In each case we varied m va[ over a wide range to cover 
both light u,d quarks as well as the heavier strange quark. Details of our simulations 
as well as the technical information on how the thermal expectation values of O were 
evaluated are in Refs. 

3. RESULTS 
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Figure 1. x/Xfft as a function of temperature for various valence quark masses. 

Based on our tests |l[ of volume dependence, made by varying N s from 8 to 16, we 
chose N s = 12. Choosing N t = 4, and making simulations at the known !3 c (N' t ) for N' t = 
6,8, 12, where f3 c is the gauge coupling at which chiral (deconfinement) transition/cross- 
over takes place for lattices with temporal extent N' t) we obtained results at T/T c =1.5, 
2 and 3 in each case. The sea quark mass in our dynamical simulations corresponds to 
14-17 MeV from the existing estimates of T c for 2-flavour QCD. Fig. [I| shows our 
results, normalized to the free field values on the same size lattice, Xfft, in the quenched 
approximation (interpolated continuous lines) and for 2 light flavours (data points) for a 
range of m va i/T c indicated on the left. Due to our choice of normalization, and the fact 
that the contribution of 2 to Xs turns out to be negligibly small for T > T c , Fig. [I] 
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shows both xs an d Xs- 

Although T c differs in the quenched and 2-ffavour QCD by a factor of 1.6-1.7, the 
respective susceptibilities in Fig. |] change by at most 5-10% for each m m z- This suggests 
that making further the strange quark dynamical may not change the results significantly. 
Using a wide range for strange quark mass of 75 to 170 MeV, the strangeness susceptibility 
can be read off from the shaded region. It shows a strong temperature dependence and is 
smaller by about 40% compared to its ideal gas value near T c . This has implications for 
phenomenology of particle abundances. Finally, the ratio x/Xfft is seen to be 0.88 (0.85) 
for 2-flavour (quenched) QCD at the smallest m va i and highest temperature we studied, 
with a mild temperature variation in its vicinity. 
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Figure 2. x/Xfft as a function of a s /n for Nf dynamical massless quarks. 



Weak coupling expansion [T(| yields x/Xfft = 1 — 2^[1 — 4y -^(1 + -g 2 -)]- Fig. |2| shows 
the predictions for various Nf along with the leading order iVj-independent prediction. 
Using a scale 2ttT for the running coupling and T c /A-jjg = 0.49(1.15) for the Nf = 2(0) 
theory ||, the values T/T c =1.5 and 3 are marked on the figure as the second (first) set. 
As one can read off from the figure, in both cases the ratio decreases with temperature in 
the range up to 3T C whereas our results in Fig. [| display an increase. Furthermore, the 
perturbative results lie significantly above in each case, being in the range 1.027-1.08 for 
2-flavour QCD and 0.986-0.994 for quenched QCD. This calls for clever resummations of 
perturbation theory or non-perturbative physics. 

One known indicator of non-perturbative physics in the plasma phase is the screening 
length in the channel with quantum numbers of pion. While it exhibits chiral symmetry 
restoration above T c by being degenerate with the corresponding scalar screening length, 
its value is much smaller than the free field value unlike that for other screening lengths. 
Fig. [3] shows X3 an d Xir (defined as a sum of the pion correlator over the entire lattice) 
as a function of M n /T. It suggests the non-perturbative physics in the two cases to be 
closely related, if not identical. 




Figure 3. 4%3/T 2 (open symbols) and Xtt/IOT 2 (filled symbols) as a function of M n /T at 
2T C (circles) and 3T C (boxes). 
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